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Abstract 

Reduced Ermakov systems are defined as Ermakov systems restricted 
' to the level surfaces of the Ermakov invariant. The condition for Lie 

. point symmetries for reduced Ermakov systems is solved yielding four 

I infinite families of systems. It is shown that SL(2, R) always is a group 

' of point symmetries for the reduced Ermakov systems. The theory is 

applied to a model example and to the equations of motion of an ion 
under a generalized Paul trap. 
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1 Introduction 

Ermakov systems play an important role in a variety of physical and 

mathematical situations. The most recent analyses involving Ermakov sys- 
tems deal with Bose-Einstein condensates and cosmological models |l]-[n], 
nonlinear supersymmetric Darboux transformations 7J, the free fall of a 
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quantum bouncing ball |HI, conformal quantum mechanics jU] and general- 
ized Hamiltonian structures ^U]- From the theoretical viewpoint, Ermakov 
systems always admit a constant of motion, the Ermakov invariant, and are 
amenable to a nonlinear superposition law ^1]. In addition, Ermakov sys- 
tems are linearizable under broad circumstances [T^ I13j . 

As is well known, the group theoretic analysis of a dynamical system is 
a subject of great relevance jHj not only for the reduction of order and the 
search for invariants for the system but also for a better understanding of its 
structural properties. The point symmetry group of Ermakov systems has 
been identified as the SL{2, R) group jlHI-Iini- More recently, using the con- 
verse to Noether's theorem, it has been shown that the Ermakov invariant 
can be associated to a dynamical symmetry, in the cases where the Ermakov 
system admits a variational formulation [20] • In addition, SL{2,R) has also 
been found |2]jj as the symmetry group of Kepler- Ermakov systems [22]) 
which can be viewed either as perturbations of the planar Kepler problem 
or of the classical Ermakov system. The purpose of this paper is to follow 
this trend from a different perspective and study the Lie point symmetries of 
Ermakov systems restricted to manifolds where the Ermakov invariant has 
a fixed constant value. The importance of this study may not be underesti- 
mated since the existence of the Ermakov invariant is automatic. 

In polar coordinates, the Ermakov system reads 

r-re^ + ^V = ^, (1) 

re + 2re = (2) 

where F and G are arbitrary functions of the angle and u, in principle, can 
depend arbitrarily on the dynamical variables. More often, is a function 
of time only, in which case it has the interpretation of a time-dependent 
frequency. For simplicity, we consider the case du/dr = 0. No mater the 
special form of a;, the Ermakov systems always possess the constant of motion 

I = \{r'9f- f G{ct>)dct>, (3) 

the so-called Ermakov invariant PP-[SI- 

The existence of a constant of motion allows the reduction of the order 
of the system. More exactly, the fourth-order system (0|2|) can be rewritten 



2 



as 

f = Air,e,t), (4) 
r^9 = B{9), (5) 

where the functions A and B are defined according to 

A{r,e,t) = -uj''r + ^{F{e) + 2I + 2 j' G{(P)d(P), (6) 

B{e) = V2il + f G{(l))d(l)f'\ (7) 

The fact that du/dr = ensures the indicated functional dependence of A 
in Notice, however, that u can freely depend, for instance, on 6, since 
this dependence can be eliminated through (jSJ. The choice du/dr = has a 
decisive infiuence on the simplification of the symmetry analysis so that we 
do not claim our results are the most general. 

Equations are a third-order dynamical system which we call the 

reduced Ermakov system. Even if we do not show explicitly, notice that the 
reduced equations do depend parametrically on I. Also, the function B is 
not identically zero except in the trivial case / = G{9) = 0, which we do not 
consider here. 

The purpose of this work is to perform the Lie point symmetry analysis 
of reduced Ermakov systems. Since for Lagrangian Ermakov systems the 
Ermakov invariant is directly related to a dynamical Noether symmetry [20] , 
it is to be expected that the algebra sl{2, R) do play a fundamental role on 
the reduced Ermakov system. 

The organization of the paper is as follows. In Section 2, the general 
symmetry conditions to be satisfied by reduced Ermakov systems and his 
generator of symmetries are determined. In Section 3, we solve the sym- 
metry conditions in the case of transformations of the time parameter not 
involving the dynamical coordinates. This yields four classes of reduced 
Ermakov systems admitting Lie point symmetries. In Section 4, we show 
that, irrespective of its specific form, the reduced Ermakov systems always 
do admit SL{2, R) as a symmetry group. This shows that, more properly, 
SL{2, R) is the group of symmetries for reduced Ermakov systems, the re- 
duction being possible thanks to a dynamical Noether symmetry in the case 
of Lagrangian Ermakov systems. In Section 5, we apply the whole formal- 
ism to a particular example, looking for reduced Ermakov systems admitting 
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a chosen Lie symmetry. In Section 6, we start from the equations for an 
ion under a generahzed Paul trap and find the circumstances under which 
these equations, written as a reduced Ermakov system, do possess Lie point 
symmetries. Section 7 is devoted to the conclusions. 

2 Lie Symmetries 

Consider infinitesimal point transformations of the form 



for functions R, S and T to be determined and infinitesimal parameter e. 
The procedure for computing Lie symmetries is well known [14j and we limit 
ourselves to sketch the critical steps in our case. The above infinitesimal 
transformation will be a Lie point symmetry of the reduced Ermakov system 
if and only if (jHEI) remains formally invariant under (jSlll(J|) up to first order 
in e, in the solution set of the reduced Ermakov system. This symmetry 
condition will imply the vanishing of two separate polynomials in r, one 
associated with the radial equation (jH), the other associated with the angular 
equation (0). Assuming the vanishing of the coefficient of different powers of 
f we get the following set of linear, coupled partial differential equations. 



r 



e 



t 



r + eR{r,9,t), 
9 + eS{r,6,t), 
t + eT{r,e,t), 



(8) 
(9) 
(10) 



(11) 



d^R 








2B d'^T ^ d'^T 2B dT 



(12) 



2B d^R 

J.2 Qj-QQ 



+ 



drdt 36 dr r"' 36 

B^ d^T 2B d^T d^T 



(13) 



^4 QQ2 ^2 QQQ^ Q^2 
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d^R BB' dR 2B d'^R d'^R 

In equations (fT^ and (fTBj) and in the sequel, a prime denotes total derivative 
with respect to 9, so that B' = dB/d9. The solutions of the equations (jllH16|) 
determine the Lie point symmetries of the reduced Ermakov system as well 
as the classes of admissible functions A and B. In the following Section we 
show four categories of solutions for the determining equations. 



3 Exact solutions 

A closer examination of ()13|) shows that the function A which specifies the 
dynamics of the radial variable in the reduced Ermakov system, will soon 
become, to some extent, determined if dT/dr ^ 0. Indeed, (fTT|) and (fT^ 
immediately give the r dependence of R and T, which, in turn, will determine 
the r dependence of A through (fT^ . Furthermore, using p4lll5|) shows that, 
for dT/dr 7^ 0, A contains only two terms, one proportional to r, the other 
to r~^. To avoid this excessively constrained situation we must have 

dT 

^-0, (17) 

a condition assumed throughout this paper. 
If (Hlj) is valid, the solution for (I11II15I) is 

R = (p(t)p(t) + r(^^))r, (18) 

S = S{e) =(yK-2fd<P ^ j 5(0) , (19) 

T = p\t), (20) 

where p is an arbitrary function of time, F is an arbitrary function of 9 and 
K a numerical constant. 

Until now, no constraint was imposed on the functions A or i? of the 
reduced Ermakov system, but there still remains the symmetry condition 
(fT^ . Inserting ()18H20|) into (fT^ . we get the following determining equation, 

UA= {-3pp + T)A + {pp + 3pp) r + ^{BT" + BT) , (21) 
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where U is the generator of Lie point symmetries for reduced Ermakov sys- 
tems, 

f-(p^ + r).| + 4 + 4. (22) 

The generator U, contains two arbitrary functions, p{t) and r(^), and one 
arbitrary numerical constant, k, imphcit in the definition (fTT^ of S{6). 
The sl{2, R) algebra is obtained from (j^^ in the particular case 

p2 = Co + cit + C2t^ , T = K = 0, (23) 

where Cq, Ci and C2 are arbitrary numerical constants. The three generators of 
the s/(2, R) algebra is obtained by taking separately each of these constants 
non zero. In comparison with the generator of point symmetries of the non 
reduced Ermakov system the new ingredients of U are r(^) and k. 
Also, p^it) is not necessarily a second degree polynomial in t. In Section 
4, we present a more detailed account of the relation between the point 
symmetries of the reduced Ermakov systems and the sl{2,R) algebra. 

Equation (j2ip can be viewed either as an equation for A or for B. We 
feel it is more productive to think of it as a determining equation for A, since 
B participates in the generator through the definition ()19|) . Following this 
choice, we find four classes of solutions for A satisfying (j^U, listed bellow. 
All these solutions are build using the differential invariants of the operator 
U, that is, the independent functions Ji and I2 for which f/(/i) = U{l2) = 0. 

3.1 The p 7^ 0, 5 7^ Case 

In this situation, the method of characteristics yields the following differential 
invariants for the generator U, 

h = -expf- r4^r(0(r;/i))l . (25) 



In fl25|) . 6 = 6(t; Ii) is a function of t as given locally by the implicit function 
theorem through (j24p . The differential invariants can be used to construct 
the solution for (j^D). The result is 

pi /■* rfr) ~ 1 /■* T(t) 

A = ^r + — exp(/ ^y^c/r) A(/i,/2) + ^exp(4 / -^dr) x 
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X 




where 6', in the integrals, is taken as a function of t through and the 
imphcit function theorem. A is an arbitrary function of the differential in- 
variants Ji and I2. 

To summarize, for any function B, the reduced Ermakov system (jUlHI) has 
a Lie point symmetry with generator (j22j) with p 7^ and 5 7^ provided 
A can be cast in the form While B remains completely arbitrary, A 

belongs to a large class of functions, including the arbitrary functions p, F 
and A and the numerical constant k. Notice that the Ermakov invariant 
enters as a parameter in the symmetry generator as well as in the function 
A. This is no surprise since the Ermakov invariant was used to eliminate 9 
from the equations of motion. 



3.2 The p ^ 0, 5 = Case 

Now the differential invariants for U are 

Ii = -, l2 = 0, (27) 
P 

and the corresponding solution for (PT|) is 

A = Pr + -^AC-,9), (28) 
p p 

where A is an arbitrary function of the indicated arguments. This class of 
solutions contains the arbitrary functions p and A, subject to p 7^ 0. Also, 
B remains completely free and does not appear in the generator. 



3.3 The = 0, r 7^ Case 

The differential invariants are 



h = \k-2 / d<P^^ r, h = t, (29) 



while 



A . ^-mL-2fi,m- r,.)..M,/.) (30) 



\ J B{(j)) 
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for A an arbitrary function depending on the differential invariants of the 
symmetry generators. Now there are the free functions F, B (which appears 
in the symmetry generator) and A, besides the numerical constant k. 



3.4 The p = r = 0, k^O Case 



This is the most simple situation. The differential invariants are r and t and 
the solution for ()21|) is A = A{r,t). The symmetry simply reflects the fact 
that ^-independent equations are invariant under rotations. 

4 Connection with the S'L(2, R) Group 

SL{2, R) is the Lie point symmetry group of non reduced Ermakov systems 
and we have to investigate the role of this transformations group 
for the reduced Ermakov systems. For simplicity, in this Section we consider 
Ermakov systems containing frequencies depending only on time. In this 
case (see equation (jH))) A is given by 



for a time-dependent frequency Q{t). Inserting (|3ip in the symmetry condi- 
tion (j2H), we get 



for S defined in (|19|) . Because (j32p has to be satisfied for arbitrary r, it can 
be split in two parts, one corresponding to r, the other to r~^, 



Equation (j34j) may be cast into a simpler form using a new independent 
variable ip given by 



A{r, e, t) = -n'{t)r + -{F{9) + B\9)) , 



(31) 



(pp + ^f)p + m'^pp + 2n(lp^)r + +{B'^T" + BB'T' + 
+ 4(F + 5^) F - {F' + 2BB') S)r-^ = , 



(32) 



pp + 3pp + A^fpp + 2VtVtp^ 
B^T" + BBT + 4(F + B^) F 



0, 

{F' + 2BB')S. 



(33) 
(34) 




(35) 
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and a new dependent variable W defined by 

„ dW 



(36) 



dip 

Using the definition (jl9|) of S as well as introducing 

H = F + B^, (37) 

equation is converted into 

+ 4H— + 2—W - K— (38) 

dif^ dip dip dip 

Equations and (jHS|) can be used to search for symmetries of specific 
reduced Ermakov systems. However, for arbitrary reduced Ermakov systems, 
that is, for completely arbitrary functions H in (jHHj) . the only possibility is 



W = K = Q. (39) 

This will be the choice if we are interested in symmetries valid for all non 
reduced Ermakov systems, regardless this specific form. Notice that W = Q 
implies F = in the symmetry generator. 

The other condition which remains is ()33|1 . which can be integrated once 
yielding Finney's [^H! equation, 

where c is a constant. However, this is a nonlinear equation, and a more 
fruitful approach for the study of symmetries is the linearizing transform 

a = p\ (41) 

giving 

a + m^a + mVL a = . (42) 

According to ()22p . the solution for this equation determines the symmetry 
generator 

d ar d 
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Now, with the rescahng 



a = a/il), T = / dfi/ip {fi), (44) 



where ip is any particular solution for the time- dependent harmonic oscillator 
equation 

^ + n'^iP = 0, (45) 

we transform (j42|) into 

The general solution is (compare with (j2Sl)) 

a = Co + ciT + C2r^ , (47) 

for constants Cq, Ci and C2- Taking separately each of these constants non zero, 
we obtain three symmetry generators for arbitrary reduced Ermakov systems. 
In the original, non rescaled variables and using ()43|) . the corresponding 
symmetry generators are 

Uo = V^'^ + ^^r— , (48) 
,2/ /■* df^ ^2 d 



Calculating the Lie brackets, the result is 

[Uo,U,]=Uo, [f/o, f/2] = 2f/i , [U,,U2] = U2, (51) 

which is the sl{2, R) algebra. This shows that the symmetry group for arbi- 
trary reduced Ermakov systems is SL{2, R). It is interesting to note that the 
algebra of the vector fields f/o, Ui and f/3 is sZ(2, R) regardless the form of 
ip (does not need to be a solution of a time-dependent harmonic oscillator). 
In addition, we notice that Uq, Ui and U2 do not depend on the Ermakov 
invariant, being generators of point transformations also in the non reduced 
space. 
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5 An Illustrative Example 



Let us consider a non reduced Ermakov system with 



J^2 gjj_|^ Q 

GiO) = 7-- , L = constant , (52) 

cos'^ 



in equation 0, producing an Ermakov invariant 

1 



and a reduced Ermakov system with 

We look for functions A{r,6,t) in the reduced Ermakov system that lead to 
Lie point symmetries as described by the class of solutions in subsection 3.1. 
Other classes of solutions may also be studied, but we restrict ourselves to 
this case in view of its generality. In fact, the solutions described in 3.1 have 
too many arbitrary functions, and we restrict ourselves to the choices 

K = 0, r(6') = Fq sin 6* , p(t) = costuot , (55) 

for constants Tq and Uq. 

In order to explicitly write the solution A, we have first to obtain the 
differential invariants (|24H25|) of the Lie symmetry. According to (|19|) and 
()54II55|) . we have 

which, substituted in fl24j) . gives the differential invariant 

/. - ^ ^ arctanh ( ^^L= \ - ^^""""^ (^7) 

2ro ^l-LV2/ \^1-Ly2l) ^0 

To obtain the second differential invariant from (j25|) . we have to solve (j57|) 
for 6, which in this case yields 



^ = arcsinh/l-^tanh 210^1 -^y(^^ + /i)l ] . (58) 
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With this result and performing the necessary quadratures, we find from (|25|) 
the differential invariant 

Again using (jSHjl and performing the necessary integrations in ()26|). we find 
the set of admissible A functions given by 

^ (1 - LV2/) sin^ / r t,, 

(1 — L'^ / 21 — sm t^j r'^ cos* Uot 

for arbitrary A depending on the differential invariants at (j57p and (j59p . We 
have redefined the arbitrary function A in ()26|) by making A A/ 12 for the 
sake of a better notation. 

To summarize, reduced Ermakov systems (^MSJ with A given by ()6Up and 
B given by do possess Lie point symmetries with generator 

d 

U = (—ujQsmuJotcosujQt + TQsm9)r— — h (61) 

or 

2ro ■ 2n^ d 2 d 

+ ^(1 - — - sm 9) — + cos ujQt — , 

cost* 2/ oO at 

as follows from ()22j) and ()55ll5fi|l . A remark applicable here and in most 
of other cases is that the Ermakov invariant does appear in the generator 
of symmetries. This is no surprise since the reduced Ermakov system was 
written after restricting the trajectories to the level surfaces of the Ermakov 
invariant. Therefore, the reduced Ermakov systems contains J as a param- 
eter. Hence, (|61|) is the generator of a point transformation on the reduced 
space and of a dynamical transformation on the non reduced space, where 
the Ermakov invariant is written as a function of r, 6 and 9 as in 



6 Application to the Generalized Paul Trap 

In practical applications, more often one does not assume the form of the 
generator and try to obtain the equations of motion for which it is a sym- 
metry, like in Section 5. Instead, one has some specific equations of motion 
and then look for symmetries. This is the approach we follow in this Sec- 
tion, searching for Lie point symmetries for the following class of Ermakov 
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systems, written in cartesian coordinates: 



where and q are initially arbitrary functions of the indicated arguments, 
and L is a constant. For constant Q and q, these are the equations of motion 
for an ion in the presence of a Paul trap |21| with equal secular frequencies. 
In this context, we call ()62H63p the equations of motion for a generalized 
Paul trap. Paul traps are a standard configuration used in ion trapping ex- 
periments 1221 • Also notice that, for q depending only on y/x, equations 
()62II63|) is a particular case of the Kepler-Ermakov systems, which are lin- 
earizable through point transformations ^3 122] • Here we ask for the classes 
of functions Q and ^ for which the corresponding reduced Ermakov systems 
do admit Lie point symmetries. 

In polar coordinates, the generalized Paul trap equations can be cast in 
the Ermakov form (PI2)) with 

f(») = 4^. = = (64) 

cos^ a cos-* a 

Notice the generalized character of a;, which is not necessarily a function of 
time only. 

The associated Ermakov invariant is 

and the reduced Ermakov system is constructed with the functions 

A(r,^,t) = -^^'Wr + ^ + ^, (66) 

B{e) = ^(l-|^-sin2^)V2, (67) 
cos a 21 

where we have defined 

a{e) = ^{B{e),tane). (68) 

Let us search for symmetries of the type shown in Section 3. For this, 
a convenient approach is to substitute A and B in the symmetry condition 
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(|21|) looking for some symmetry generator. The symmetry condition (|2H) 
then gives 

{pp + 3pp + An^pp + 2nnp'^)r + {-S{e)a' + {3T{e) - pp)a)r-'^ (69) 
,2/(1 -LV2/-sin2^) ^„ L^sin^^, , 

COS"^ U COS'^ u 

where S{9) depends on T{9) according to (fTUj) . Compare with (jSH). 
There are similarities since B{6) here is the same as in Section 5, but now 
there is a contribution proportional to r~^, peculiar to Kepler- Ermakov sys- 
tems. 

As in Section 5, we split the symmetry condition ()(i9j] into three equations, 
corresponding to different powers of r, 

pp + 3pp + 4Q^pp + 2QQp^ = 0, (70) 
Sie)a'+{pp-3T{e))a = 0, (71) 
2/(1 - LV2/ - sin^ 9) ^„ sin 9 ^, 

— '-^ ^r" - — — r' + 8/r = o. 72 

cos^ 9 cos'^ 9 

For consistency with equation ()68p . in equation ()71|) we must have 

PP = ^. (73) 
for a constant Qq. Integrating, we get 

p = (pI + 2notY/\ (74) 



where po is a constant. Inserting this into (fTOj) . we obtain the following class 
of frequencies Q, 

pi + 2not 

Equation (f7^ can be best handled using the new independent variable 



dp 1 sin^ 
if = / ^, , = -^= arcsmf , , (76) 

' B{fi) V2I VI-LV2/ 



yielding 



+ 8/r = o, (77) 
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with solution 

r = Ti cos(2 V27(^) + T2 sm{2V2Iip) , (78) 

where Fi and r2 are constants. 

There remains the equation (f7T|). Calculating S{9) via (fT^ and using the 
new independent variable ip, this equation reads 

(«:+^)^ + (fio-3r)a = 0, (79) 

with r as in (fTHj) . 

There is no difficulty at solving (f?^ in the general case, but the solution 
is somewhat complicated. It seems more illustrative to obtain the general 
solution in some specific cases, listed bellow. 

6.1 The ^ 0, Ti = r2 = Case 

The general solution for (f7^ is 

f ^0 ■ ( sine 
(T = (Toe « =(Joexp( =arcsm( ^=)) , (80) 



for any constant (Tq. 



6.2 The rf + ri 7^ 0, «: = Qo = Case 

In this situation the general solution is 

where (Tq is a constant and F is given by (fTHj). with (y9 defined in (fTBjl . 

In both subcases ()80|) and (jHT|) . the solution contains the Ermakov invari- 
ant which is dependent on Hence, at least in the generalized Paul 
trap case, we have not found solutions for which o" is a function of 9 only, 
characterizing a Kepler-Ermakov system. 

To summarize, we found the functions a (alternatively, ^) and VL so that 
the reduced Ermakov system for the generalized Paul trap equations do pos- 
sess Lie point symmetries. The functions a do satisfy (f?^ . as in subsections 
6.1 and 6.2, while the frequency VL is of the form (f73j) . The corresponding 
generator of symmetry is build with p given by (j74j) and F given by (j78j) . 
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7 Conclusion 



We presented the general treatment for Lie point symmetries of reduced 
Ermakov systems. We found four classes of reduced Ermakov systems pos- 
sessing Lie point symmetries, all of them involving arbitrary functions. We 
have applied the results to two different systems, including a generalized 
Paul trap. From the theoretical viewpoint, the most important result we 
have found is the fact that the SL{2, R) group is more exactly a property 
of the reduced Ermakov system. For Lagrangian Ermakov systems, the exis- 
tence of the Ermakov invariant follows from a dynamical symmetry. Then, 
for reduced Lagrangian Ermakov systems, the symmetry structure can be 
separated in two distinct parts: a dynamical symmetry leading to the Er- 
makov invariant, and the SL{2, R) group for the reduced Ermakov system, 
the later a consequence of the first symmetry. 
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